CERN-PH-TH/2010-004, UMN-TH-2831/10, FTPI-MINN-10/02 



Supersymmetric Dark Matter Candidates 

John Ellis a and Keith A. Olive b 

a TH Division, Physics Department, CERN, Geneva, Switzerland 
b William I. Fine Theoretical Physics Institute, University of Minnesota, 
Minneapolis, MN 55455, USA 

Abstract 

After reviewing the theoretical, phenomenological and experimental motiva- 
tions for supersymmetric extensions of the Standard Model, we recall that su- 
persymmetric relics from the Big Bang are expected in models that conserve R 
parity. We then discuss possible supersymmetric dark matter candidates, focusing 
on the lightest neutralino and the gravitino. In the latter case, the next-to-lightest 
supersymmetric particle is expected to be long-lived, and possible candidates in- 
clude spartners of the tau lepton, top quark and neutrino. We then discuss the 
roles of the renormalization-group equations and electroweak symmetry breaking 
in delimiting the supersymmetric parameter space. We discuss in particular the 
constrained minimal extension of the Standard Model (CMSSM), in which the 
super symmetry-breaking parameters are assumed to be universal at the grand 
unification scale, presenting predictions from a frequentist analysis of its param- 
eter space. We also discuss astrophysical and cosmological constraints on grav- 
itino dark matter models, as well as the parameter space of minimal supergravity 
(mSUGRA) models in which there are extra relations between the trilinear and bi- 
linear supersymmetry-breaking parameters, and between the gravitino and scalar 
masses. Finally, we discuss models with non-universal supersymmetry-breaking 
contributions to Higgs masses, and models in which the supersymmetry-breaking 
parameters are universal at some scale below that of grand unification. 
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Gianfranco Bertone Copyright 2010 Cambridge University Press. Chapter 8, pp. 
142-163 Hardback ISBN 9780521763684, 
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1 Motivations 



Supersymmetry is one of the best-motivated proposals for physics beyond the 
Standard Model. There are many idealistic motivations for believing in su- 
persymmetry, such as its intrinsic elegance, its ability to link matter particles 
and force carriers, its ability to link gravity to the other fundamental interac- 
tions, its essential role in string theory, etc. However, none of these aesthetic 
motivations gives any hint as to the energy scale at which supersymmetry 
might appear. The following are the principal utilitarian reasons to think 
that supersymmetry might appear at some energy accessible to forthcoming 
experiments. 

The first and primary of these was the observation that supersymmetry 
could help stabilize the mass scale of electroweak symmetry breaking, by 
cancelling the quadratic divergences in the radiative corrections to the mass- 
squared of the Higgs boson [TJ |2j [3], and by extension to the masses of other 
Standard Model particles. This motivation suggests that sparticles weigh 
less than about 1 TeV, but the exact mass scale depends on the amount of 
fine-tuning that one is prepared to tolerate. 

Historically, the second motivation for low-scale supersymmetry, and the 
one that interests us most here, was the observation that the lightest super- 
symmetric particle (LSP) in models with conserved R parity, being heavy and 
naturally neutral and stable, would be an excellent candidate for dark mat- 
ter jHEj. This motivation requires that the lightest supersymmetric particle 
should weigh less than about 1 TeV, if it had once been in thermal equilib- 
rium in the early Universe. This would have been the case for a neutralino \ 
or a sneutrino v LSP, and the argument can be extended to a gravitino LSP 
because it may be produced in the decays of heavier, equilibrated sparticles. 

The third reason that emerged for thinking that supersymmetry may be 
accessible to experiment was the observation that including sparticles in the 
renormalization-group equations (RGEs) for the gauge couplings of the Stan- 
dard Model would permit them to unify [6j [7J El [9J [TO], whereas unification 
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would not occur if only the Standard Model particles were included in the 
RGEs. However, this argument does not constrain the supersymmetric mass 
scale very precisely: scales up to about 10 TeV or perhaps more could be 
compatible with grand unification. 

The fourth motivation is the fact that the Higgs boson is (presumably) 
relatively light, according to the precision electroweak data - an argument 
reinforced by the negative results (so far) of searches for the Higgs boson at 
the Fermilab Tevatron collider. It has been known for some 20 years that 
the lightest supersymmetric Higgs boson should weigh no more than about 
140 GeV, at least in simple models PU 02 H31 H H5] ■ Since the early 1990s, 
the precision electroweak noose has been tightening, and the best indication 
now (incorporating the negative results of searches at LEP and the Tevatron) 
is that the Higgs boson probably weighs less than about 140 GeV [TU ITT] , 
in perfect agreement with the supersymmetric prediction. 

Fifthly, if the Higgs boson is indeed so light, the present electroweak 
vacuum would be destabilized by radiative corrections due to the top quark, 
unless the Standard Model is supplemented by additional scalar particles [T8] . 
This would be automatic in supersymmetry, and one can extend the argu- 
ment to 'prove' that any mechanism to stabilize the electroweak vacuum must 
look very much like supersymmetry. 

There is a sixth argument that is still controversial, namely the anomalous 
magnetic moment of the muon, g> M — 2. As is well known, the experimental 
measurement of this quantity [19] disagrees with the Standard Model pre- 
diction |20j, if this is calculated using low-energy e + e~ annihilation data. 
On the other hand, the discrepancy with the Standard Model is greatly re- 
duced if one uses r decay data to estimate the Standard Model contribution 
to gfi — 2. Normally, one would prefer to use e + e~ data, since they are re- 
lated more directly to — 2, with no need to worry about isospin violations, 
etc. Measurements by the BABAR collaboration using the radiative-return 
method [21] yield a result intermediate between the previous e + e~ data and 
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t decay data. Until the discrepancy between these data sets have been ironed 
out, one should take g M — 2 cum grano salis. 

2 The MSSM and R parity 

We refer to [22j E3] for the general structure of supersymmetric theories. 
We restrict ourselves here to theories with a single supersymmetry charge, 
called simple or N = 1 supersymmetry, as these are the only ones able to 
accommodate chiral fermions and hence the violation of parity and charge 
conjugation. We recall that the basic building building blocks of N = 1 su- 
persymmetric models are so-called chiral supermultiplets, each consisting of a 
Weyl fermion and a complex scalar, and gauge supermultiplets, each consist- 
ing of a gauge field and a gaugino fermion. The renormalizable interactions 
between the chiral supermultiplets are characterized by a superpotential that 
couples the chiral supermultiplets in bilinear and trilinear combinations that 
yield masses and Yukawa interactions, and by gauge interactions. In this 
framework, bosons and fermions must appear in pairs with identical internal 
quantum numbers. Since the known particles do not pair up in this way, 
it is necessary to postulate unseen particles to partner those known in the 
Standard Model. 

In order to construct the minimal supersymmetric extension of the Stan- 
dard Model (MSSM) [2H [251 12H], one starts with the complete set of chiral 
fermions needed in the Standard Model, and adds a complex scalar super- 
partner to each Weyl fermion, so that each matter field in the Standard Model 
is extended to a chiral supermultiplet. These are denoted by L\ Q\ e c , d c and 
u c , where i,j are SU(2)i doublet indices and generation indices have been 
suppressed as were color indices for the quarks. In order to avoid a triangle 
anomaly, Higgs supermultiplets must appear in pairs with opposite hyper- 
charges, and the minimal possibility is a single pair H^H^. One must also 
add a gaugino for each of the gauge bosons in the Standard Model so as to 
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complete the gauge supermultiplets. The minimal supersymmetric standard 
model (MSSM) [27J is defined by this minimal field content and the minimal 
superpotential necessary to account for the necessary Yukawa couplings and 
mass terms, namely: 

W = e l3 {y e H{L l e c + y d H 3 1 Q i d c + y u E\Q'u c ) + e^tf^. (1) 

In (flj) , the Yukawa couplings, y, are all 3 x 3 matrices in generation space, 
with no generation indices for the Higgs multiplets. A second reason for 
requiring two Higgs doublets in the MSSM is that the superpotential must 
be a holomorphic function of the chiral superfields. This implies that there 
would be no way to account for all of the Yukawa terms for both up- and 
down- type quarks, as well as charged leptons, with a single Higgs doublet. 
The physical Higgs spectrum then contains five states: two charged Higgs 
bosons H ± , two scalar neutral Higgs bosons h,H, and a pseudoscalar Higgs 
boson A. The final bilinear mixing term in ([1]) must be included in the 
superpotential, in order to avoid a massless Higgs state. 

The MSSM must be coupled to gravity, which requires the introduction 
of a graviton supermultiplet containing a spin-3/2 gravitino as well as the 
spin-2 graviton itself, which may or not be coupled minimally to the MSSM. 
The consistency of supergravity at the quantum level requires the breaking of 
supersymmetry to be spontaneous, with the gravitino mass acting as an or- 
der parameter [281 12H] . The mechanism whereby supersymmetry is broken is 
unknown, as is how this feeds into the MSSM. We adopt here a phenomeno- 
logical approach, parametrizing the results of this mechanism in terms of 
differing amounts of explicit supersymmetry breaking in the masses and cou- 
plings of the unseen supersymmetric partners of Standard Model particles 
[301I3HE3]. 

In order to preserve the hierarchy between the electroweak and GUT or 
Planck scales, it is necessary that this explicit breaking of supersymmetry 
be 'soft', i.e., in such a way that the theory remains free of quadratic diver- 
gences, which is possible with the insertion of weak scale mass terms in the 
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Lagrangian (32]- The possible forms for such terms are 



Csoft = - l -M a \ a \ a -I( m 2 )«<^* 



(2) 



where the M a are masses for the gauginos A a , m 2 is a matrix of soft scalar 
masses-squared that carries two field indices, a, (3, for scalars <p a , A is a trilin- 
ear coupling term with three field indices, and B is a bilinear supersymmetry 
breaking term associated with a superpotential bilinear mass term such as 
in Eq. [U Masses for the gauge bosons are, as usual, induced by the spon- 
taneous breaking of gauge invariance, and the masses for chiral fermions are 
induced by the Yukawa superpotential terms when the electroweak gauge 
symmetry is broken. For a more complete discussion of supersymmetry and 
the construction of the MSSM see [33l El ESI |36] . 

In defining the MSSM, we have limited the model to contain a minimal 
field content: the only new fields are those which are required by supersym- 
metry Consequently, apart from superpartners, only the Higgs sector was 
enlarged from one doublet to two. Moreover, in writing the superpotential 
([1]), we have also made a minimal choice regarding interactions. We have 
limited the types of interactions to include only the minimal set required in 
the Standard Model and its supersymmetric generalization. 

However, even with the minimal field content, there are several other 
superpotential terms that one could envision adding to ([I]) which are consis- 
tent with all of the gauge symmetries of the theory. Specifically, one could 
consider adding any or all of the following terms that violate i?-parity: 



Each of the terms in ((3]) has one or more suppressed generation indices. We 
note that the terms proportional to A, A', and ^! both violate lepton number 
by one unit, whereas the term proportional to A" violates baryon number by 
one unit. 



W R 




(3) 
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Each of the terms in ([3]) predicts new particle interactions and can be 
to some extent constrained by the lack of observed exotic phenomena. In 
particular, any combination of terms which violate both baryon and lepton 
number would be unacceptable, unless the product of coefficients was ex- 
tremely small. For example, consider the possibility that both A' and A" 
were non-zero. This would lead to the following proton decay processes: 
p —> e + 7r°, /i + 7T°, utt + , b>K + , etc. The rate of proton decay due to this process 
would have no suppression by any superheavy masses, since there is no GUT- 
or Planck-scale physics involved: this is a purely (supersymmetric) Standard 
Model interaction involving only the electroweak scale. The (inverse) rate 
can be easily estimated to be 

r; 1 ~ % ~ 10 8 GeV-\ (4) 

assuming a supersymmetry breaking scale of m of order 100 GeV. This should 
be compared with current limits to the proton life-time of > 10 63 GeV -1 . 
Clearly the product of A' and A" must be very small, if not exactly zero. 

It is possible to eliminate the unwanted superpotential terms by impos- 
ing a discrete symmetry on the theory called i?-parity [37]. This can be 
represented as 

R = (_1)3B+L+2 S; ( 5) 

where B, L, and s are the baryon number, lepton number, and spin respec- 
tively. It is easy to see that, with the definition (J5J), all the known Standard 
Model particles have i?-parity +1. For example, the electron has B = 0, 
L = —1, and s = 1/2, and the photon has B = L = and s = 1, so in 
both cases R — 1. Similarly, it is clear that all superpartners of the known 
Standard model particles have R = — 1, since they must have the same value 
of B and L as their conventional partners, but differ by 1/2 unit of spin. If 
i?-parity is exactly conserved, then all four superpotential terms in ([3]) must 
be absent. 
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The additive conservation of the quantum numbers B, L, and s implies 
that .R-parity must be conserved multiplicatively. A first important corollary 
is that the collisions of conventional particles must always produce supersym- 
metric particles in pairs, and a second corollary is that heavier supersymmet- 
ric particles can decay only into lighter supersymmetric particles. For our 
purposes here, an even more important corollary of i?-parity conservation 
is the prediction that the lightest supersymmetric particle (LSP) must be 
stable, because it has no legal decay mode. In much the same way that 
baryon number conservation predicts proton stability, i?-parity predicts that 
the lightest R — — 1 state is stable. This makes supersymmetry an extremely 
interesting theory from the astrophysical point of view, as the LSP naturally 
becomes a viable dark matter candidate (31 E] . 

3 Possible Supersymmetric Dark Matter Can- 
didates 

What options are available in the MSSM for the stable LSP? Any electrically- 
charged LSP would bind to conventional matter, and be detectable as an 
anomalous heavy nucleus, since the 'Bohr radius' for the LSP 'atom' would 
be less than the nuclear radius. Similarly, strongly-interacting LSPs would 
also form anomalous heavy nuclei. However, experiments searching for such 
objects [38j [39j HO] have excluded their presence on Earth down to an abun- 
dance far lower than the expected abundance for the LSP (see below for 
more details how this is calculated). Therefore, the stable LSP is presum- 
ably electrically-neutral and can have only weak interactions. For this reason, 
the commonly-expected signature of supersymmetric particle production at 
colliders is missing energy carried away by undetected LSPs. 

This still leaves us with several possible dark matter candidates in the 
MSSM, specifically the sneutrino with spin zero, the neutralino with spin 
1/2, and the gravitino with spin 3/2. However, a sneutrino LSP would 
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have relatively large coherent interactions with heavy nuclei, and experi- 
ments searching directly for the scattering of massive dark matter particles 
on nuclei exclude a stable sneutrino weighing between a few GeV and several 
TeV [41 J. The possible loophole of a very light sneutrino was excluded by 
measurements of the invisible Z-boson decay rate at LEP [42] . 

The LSP candidate that is considered most often is the lightest neutralino. 
In the MSSM there are four neutralinos, each of which is a linear combination 
of the following R — — 1 neutral fermions [5]: the neutral wino W 3 (the 
partner of the third component of the SU(2) L triplet of weak gauge bosons); 
the U(l) bino B\ and two neutral Higgsinos Hi and H 2 (the supersymmetric 
partners of the neutral components of the two Higgs doublets). 

The composition of the LSP x can De expressed as a linear combination 
of these fields: 

X = aB + (5W 3 + nfH x + 5H 2 , (6) 

whose mass and composition are determined by the SU(2)l and U(l) gaugino 
masses, the Higgs mixing parameter /x, and tan/3, the ratio of the 

vacuum expectation values t>i j2 =< 0|ifi )2 |0 > of the two neutral Higgs fields 
tan j3 = v%jv\. The mass of the LSP x an d the mixing coefficients a, /3, 7 and 
5 in ([H]) for the neutralino components that compose the LSP can be found 
by diagonalizing the mass matrix 
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In different regions of the supersymmetric parameter space, the LSP may be 
more bino-like, wino-like, or Higgsino-like, depending on the relative magni- 
tudes of the coefficients a, {3, 7 and S. 

The relic abundance of an LSP candidate such as the lightest neutralino 
is calculated by solving the Boltzmann equation for the LSP number density 
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in an expanding Universe: 



dn 



3-n- (av)(n 2 -n 2 ), 



(8) 



where n is the equilibrium number density of neutralinos. Defining the 
quantity / = n/T 3 , we can rewrite this equation in terms of the reduced 
temperature x = T /m x : 



where Gn is Newton's constant and N is the number of relativistic degrees of 
freedom at a given temperature. The solution to this equation at late times 
and low temperatures, and hence small x, yields a constant value of /, so 
that n oc T 3 . 

The technique [13] used to determine the neutralino relic density is sim- 
ilar to that used previously for computing the relic abundance of massive 
neutrinos [^H H5j 26] , with the substitution of the appropriate annihilation 
cross section. This and hence the relic density depend on additional param- 
eters in the MSSM beyond Mi, M 2 , /i, and tan fl, which include the sfermion 
masses, rrij and mass of the pseudoscalar Higgs boson, m^. In much of the 
parameter space of interest, the LSP is a bino and the annihilation proceeds 
mainly through crossed t-channel sfermion exchange. The exception is if the 
sum of two neutralino masses happens to lie near a direct-channel pole, such 
as m x ~ mz/2 or in which case there are large contributions to the an- 

nihilation through direct s-channel resonance exchange. Since the neutralino 
is a Majorana fermion, away from such a resonance the s-wave part of the 
annihilation cross section is generally suppressed by the outgoing fermion 
masses, and the annihilation occurs mainly through the p wave, which is also 
suppressed because the annihilating LSPs are non-relativistic at low temper- 
atures (small x). This means that one can approximate the annihilation cross 
section including p-wave corrections by incorporating a term proportional to 
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the temperature if neutralinos are in thermal equilibrium: av = a + bx + . . . , 
where the expansion coefficients a, b are model-dependent. 

Annihilations in the early Universe continue until the annihilation rate 
T ~ avn x drops below the expansion rate, after which it is a good first 
approximation to assume that annihilations are negligible - the freeze-out 
approximation. The final neutralino relic density, expressed as a fraction fl x 
of the critical energy density and denoting the present-day Hubble expansion 
rate as h in units of 100 km/s/Mpc, can be written as [5] 



where (T x /T 7 ) 3 accounts for the subsequent reheating of the photon tem- 
perature with respect to x> due to the annihilations of particles with mass 
m < Xfm x [U] HH|, and Xf — Tf/m x is proportional to the freeze-out tem- 
perature Tf. Eq. ([10]) yields a very good approximation to the relic density 
except near direct s-channel annihilation poles or thresholds, and in regions 
where the LSP is nearly degenerate with the next lightest supersymmetric 
particle |49j. 

When there are several particle species i that are nearly degenerate in 
mass, coannihilations between the different species become important. In 
this case [19], the rate equation ([8]) still applies, provided n is interpreted as 
the total number density, 




(10) 



n = 



(11) 



n is interpreted as the total equilibrium number density, 



(12) 



and the effective annihilation cross section as 




(13) 
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In eq. (jHJ), m x is now understood to be the mass of the lightest sparticle 
under consideration. 

We turn finally to the third LSP candidate within the MSSM, namely 
the gravitino. Since it has only gravitational-strength interactions, it is not 
expected to have been in thermal equilibrium in the early Universe. However, 
it could have been produced in high-energy particle collisions in the early 
Universe, or in the decays of heavier supersymmetric particles. The fact 
that the gravitino has only gravitational-strength interactions implies that 
only decays of the next-to-lightest supersymmetric particle (NLSP) would 
be significant sources of gravitinos, and the NLSP would be metastable. 
As we discuss in more detail later, there are important cosmological and 
astrophysical constraints on the possible mass and lifetime of the NLSP, 
derived principally from the agreement between astrophysical observations 
and Big-Bang Nucleosynthesis calculations of light-element abundances. 

What might be the nature of the NLSP be in such a gravitino LSP sce- 
nario? One option is the lighter of the two supersymmetric partners of the 
t lepton, denoted by T\. Being a metastable charged particle, it would have 
a distinctive experimental signature at the LHC or other colliders. Studies 
within such a scenario have shown that the mass of the f\ could be measured 
very accurately, and that one could easily reconstruct heavier sparticles that 
decay into the f\ [50] . 

Alternatively, the NLSP might be the lighter supersymmetric partner 
of the top quark, denoted by t\ [51] [52], [53], which would have even more 
distinctive signatures at the LHC. Immediately after production, it would 
become confined inside a charged or neutral hadron. As it moves through 
an LHC detector, it would have a high probability of changing its charge 
as it interacts with the material in the detector. This combined with its 
non-relativistic velocity would provide a truly distinctive signature. 

Yet another possibility is that the NLSP might be some flavour of sneu- 
trino [54], in which case the characteristic signature would be missing energy 
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carried away by the metastable sneutrino. This could nevertheless be dis- 
tinguished from the conventional missing-energy signature of a neutralino 
LSP (or NLSP), because the final states would be more likely to include the 
charged lepton with the same flavour as the sneutrino NLSP, either e, \i or 
r. 

These are just a few examples of the possible alternatives to the conven- 
tional missing-energy signature of supersymmetry. Studies have shown that 
the LHC would also have good prospects for detecting such signatures. 



4 Renormalizat ion- Group Equations and Elec- 
troweak Symmetry Breaking 

The fact that measurements of the strengths of the Standard Model gauge 
interactions measured at low energies are in excellent agreement with the 
predictions of a supersymmetric gauge theory 0, [7J EJ EJ [TO] was already cited 
as an important motivation for low-energy supersymmetry. It can also be 
regarded as a motivation for thinking that other parameters of the effective 
low-energy theory, e.g., the soft supersymmetry breaking parameters can 
also be calculated and related using renormalizat ion-group equations (RGEs) 
below the grand unification scale. For example, the one-loop RGEs for the 
gaugino masses are: 

^ = -ka^/An (14) 
at 

If the gaugino masses have a common value mu^ at the grand unification 
scale, these equations can be used to relate the physical low-energy, on-shell 
values of the gaugino masses to the corresponding gauge coupling strengths 

Oti. 

Miit) = *® m 1/2 , (15) 
ai(M GUT ) 



which implies that 



Mi M 2 M 3 

— = — = — ( 16 ) 
9i 92 9s 
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at the one-loop level. When applying this relation within a specific grand 
unified theory, one must remember to incorporate the difference of the nor- 
malization of the U(l) factor from that in the Standard Model, so that we 
have Mi = |^-M 2 for the one- loop relation between the bino and wino 
masses. Also, the simple relations ( TIT)]) are modified by threshold correc- 
tions at the electroweak scale, and by two-loop effects in the RGEs. The 
soft supersymmetry-breaking scalar masses- squared m 2 and the trilinear cou- 
plings A are renormalized analogously to (114p . with the difference that Yukawa 
interactions contribute as well as gauge interactions. However, the Yukawa 
contributions are small, except for the supersymmetric partners of third- 
generation fermions. 

As described above, the MSSM has over 100 undetermined parameters, 
which are mainly associated with the breaking of supersymmetry. It is often 
assumed that the soft supersymmetry-breaking parameters M a , m 2 and A 
have some universality properties. There are phenomenological arguments, 
based on the success of the Standard Model in describing the observed sup- 
pression of flavour- changing interactions, that, at some input scale (often 
assumed to be that of grand unification), the parameters m 2 and A must be 
universal for supersymmetric particles with the same gauge quantum num- 
bers, e.g., the supersymmetric partners of the e, fi and r. There is no strong 
argument why these parameters should be universal for supersymmetric par- 
ticles with different quantum numbers, e.g., d, u and e, though this may occur 
in some grand unified theories, as may unification of the gaugino masses M a . 
The simplified version of the MSSM in which universality at the grand uni- 
fication scale is assumed for each of M a , m 2 and A is called the constrained 
MSSM (CMSSM)[53[5g[53[5El[59j[6a 

uncatsiinitsitsiirzits]- 

Once one has chosen a set of boundary conditions at the grand unifica- 
tion scale and run the RGEs down to the electroweak scale, one must check 
the properties of the electroweak vacuum, which are characterized by speci- 
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fying the mass of the Z boson, Mz, and the ratio of the two Higgs vacuum 
expectation values, tan (3. These electroweak symmetry-breaking conditions 
should be used as consistency conditions on the solutions to the RGEs, e..g., 
of the CMSSM. They are frequently used to fix, as functions of the input 
values of the common gaugino mass m^, m, A and tan (3, the magnitudes of 
the Higgs mixing mass parameter, /i, and of the bilinear coupling, B, which 
determines the pseudoscalar Higgs mass, tha- The sign of /i remains free. 

An example of the running of the mass parameters in the CMSSM as 
functions of the renormalization scale is shown in Fig. (TJ using as inputs the 
choices m 1/2 = 250 GeV, m = 100 GeV, tan/3 = 3, A = 0, and /x < 0. We 
notice in the figure several characteristic features of the sparticle spectrum. 
For example, the colored sparticles are typically the heaviest, because of the 
large positive corrections to their masses arising from a3-dependent terms 
in the RGEs. Also, one finds that the bino, B, is typically the lightest 
sparticle. Most importantly, we notice that one of the Higgs masses squared, 
goes negative, triggering electroweak symmetry breaking [79j |80j Ell E2j 183] . 
(The negative sign in the figure refers to the sign of the mass squared, even 
though it is the mass of the sparticles which is depicted.) 

5 The CMSSM 

For given values of tan /3, A Q , and sgn(fi), the regions of the CMSSM param- 
eter space that yield an acceptable relic density and satisfy the other phe- 
nomenological constraints may conveniently be displayed in the {m 1 / 2 ,mo) 
plane. Fig. [2] displays, for tan/3 = 10 (a) and 50 (b), the impacts of the most 
relevant constraints. These include the LEP lower limits on the chargino 
mass: m x ± > 104 GeV [53] . on the selectron mass: mg > 99 GeV [85] and 
on the Higgs mass: rrih > 114 GeV [861 ISTj. The former two constrain 7711/2 
and m directly via the sparticle masses, and the latter indirectly via the 
sensitivity of radiative corrections to the Higgs mass to the sparticle masses, 
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Figure 1: The renormalization- group evolution of the mass parameters in the 
CMSSM, assuming m 1/2 = 250 GeV, m = 100 GeV, tan/3 = 3, A = 0, 
and fi < 0. We thank Toby Falk for providing this figure. 

principally m^. Here the code FeynHiggs [88, 89J is used for the calculation 
of nth- It would be prudent to assign an uncertainty of 3 GeV to this calcula- 
tion. Nevertheless, the Higgs limit imposes important constraints, principally 
on m.i/2 and particularly at low tan f3. Another constraint is the requirement 
that the branching ratio for b — )■ s'j be consistent with the experimental 
measurements [90J. These measurements agree with the Standard Model, 
and therefore provide bounds on MSSM particles [91], such as the chargino 
and charged Higgs bosons, in particular. Typically, the b — > 57 constraint is 
more important for \i < 0, but it is also relevant for /x > 0, particularly when 
tan /3 is large. The constraint imposed by measurements of b — > 57 also ex- 
clude small values of mm. Finally, there are regions of the (m^, m ) plane 
that are favoured by the Brookhaven National Laboratory measurement [19] 
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of — Here we assume the Standard Model calculation [2D] of g M — 2 using 
e + e~ data, and indicate by dashed and solid lines the contours of 1- and 2-er 
level deviations induced by supersymmetry. 




m 1/2 (GeV) m 1/2 (GeV) 



Figure 2: The (rni/2, Wo) planes for (a) tan/3 = 10 and tan = 50, 
assuming n > 0, Aq = 0, m t = 175 GeV and mb(mb)^ = 4.25 GeV. 
The near-vertical (red) dot-dashed lines are the contours for rah = 114 GeV, 
and the near-vertical (black) dashed line is the contour m x ± = 104 GeV. 
Also shown by the dot-dashed curve in the lower left is the region excluded 
by the LEP bound mg > 99 GeV. The medium (dark green) shaded region is 
excluded byb^- s'y, and the light (turquoise) shaded area is the cosmologically 
preferred region. In the dark (brick red) shaded region, the LSP is the charged 
f\ . The region allowed by the E821 measurement of a M at the 2-a level, is 
shaded (pink) and bounded by solid black lines, with dashed lines indicating 
the 1-a ranges. 

The most precise constraint on supersymmetry may be that provided 
by the density of cold dark matter, as determined from astrophysical and 
cosmological measurements by WMAP and other experiments |92j : 

n CDM = 0.1099 ±0.0062. (17) 

Applied straightforwardly to the relic LSP density Qisph 2 , this would give 
a very tight relation between supersymmetric model parameters, fixing some 
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combination of them at the % level, which would essentially reduce the di- 
mensionality of the supersymmetric parameter space by one unit. Let us 
assume for now that the LSP is the lightest neutralino x> whose density is 
usually thought to be fixed by freeze-out from thermal equilibrium in the 
early Universe, as discussed previously. In this case, respecting the con- 
straint (TTTj) would force the CMSSM into one of the narrow WMAP 'strips' 
in planar projections of the parameters |74J . as illustrated by the narrow light 
(turquoise) regions in Fig. [2J However, caution should be exercised before 
jumping to this conclusion. 

Supersymmetry might not be the only contribution to the cold dark mat- 
ter, in which case ( |17|) should be interpreted as an upper limit on VL^sph 2 . 
However, most of the supersymmetric parameter space in simple models gives 
a supersymmetric relic density that exceeds the WMAP range (ITTj) . e.g., 
above the WMAP 'strip' in Fig. [2j and the regions with lower density gener- 
ally correspond to /overvalues of the sparticle masses, i.e., below the WMAP 
'strip' in Fig. 

However, even if one takes them seriously, the locations of these WMAP 
'strips' do vary significantly with the choices of other supersymmetric param- 
eters, as can be seen by comparing the cases of tan/3 = 10, 50 in Fig. |2]^a, b). 
As one varies tan/3, the WMAP 'strips' cover much of the (mi/2,mo) plane. 

Several different regions of the WMAP 'strips' in the CMSSM (mi/2,mo) 
plane can be distinguished, in which different dynamical processes are dom- 
inant. At low values of mi/2 and mo, simple \ — % annihilations via crossed- 
channel sfermion exchange are dominant, but this 'bulk' region is now largely 
excluded by the LEP lower limit on the Higgs mass, m^. At larger m^, but 
relatively small mo, close to the boundary of the region where the lighter 
stau is lighter than the lightest neutralino: < m x , coannihilation be- 
tween the x an d sleptons are important in suppressing the relic x density 
into the WMAP range (JT7j) . as seen in Fig. |2j At larger mi/2, mo and tan/3, 
the relic x density may be reduced by rapid annihilation through direct- 
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channel H, A Higgs bosons, as seen in Fig. [2(b). Finally, the relic density 
can again be brought down into the WMAP range (ITT)) at large mo (not 
shown in Fig. [2]), in the 'focus-point' region close the boundary where elec- 
troweak symmetry breaking ceases to be possible and the lightest neutralino 
X acquires a significant higgsino component [53] . 

As seen in Fig. [2j the relic density constraint is compatible with rel- 
atively large values of m\/2 and itlq, and it is interesting to look for any 
indication where the supersymmetric mass scale might lie within this range, 
using the available phenomenological and cosmological constraints. A global 
likelihood analysis enables one to pin down the available parameter space 
in the CMSSM and the related models discussed later. One can avoid the 
dependence on priors by performing a pure likelihood analysis as in [HI] , or a 
purely x 2 -based fit as done in [951 EH] . Here we present results from one such 
analysis [HTJ, which used a Markov- Chain Monte Carlo (MCMC) technique 
to explore efficiently the likelihood function in the parameter space of the 
CMSSM. A full list of the observables and the values assumed for them in 
this global analysis are given in [96], as updated in [9T] . 

The 68% and 95% confidence-level (C.L.) regions in the (mi/ 2 ,mo) plane 
of the CMSSM are shown in Fig. [3] [97] . Also shown for comparison are the 
physics reaches of ATLAS and CMS with 1/fb of integrated luminosity [981 
[99] . (MET stands for missing transverse energy, SS stands for same-sign 
dilepton pairs, and the sensitivity for finding the lightest Higgs boson in 
cascade decays of supersymmetric particles is calculated for 2/fb of data.) 
The likelihood analysis assumed fi > 0, as motivated by the sign of the 
apparent discrepancy in — 2, but sampled all values of tan/3 and A : the 
experimental sensitivities were estimated assuming tan (3 = 10 and Aq = 0, 
but are probably not very sensitive to these assumptions. The global maxima 
of the likelihood function (indicated by the black dot) is at mi/2 = 310 GeV, 
m = 60 GeV, A = 240 GeV, tan/3 = 11 and x 2 /N dof = 20.4/19 (37% 
probability). It is encouraging that the best-fit points lie well within the LHC 
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discovery range, as do the 68% and most of the 95% C.L. regions. It is also 
encouraging that the two best-fit points have similar values of m^jfflo and 
tan /3, the most important parameters for the sparticle spectrum, indicating 
that the likelihood analysis is relatively insensitive to the theoretical model 
assumptions. 
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Figure 3: The (mo,™^) plane in the CMSSM showing the regions favoured 
in a likelihood analysis at the 68% (blue) and 95% (red) confidence levels [91]. 
The best-fit point is shown as the black point. Also shown are the discovery 
contours in different channels for the LHC with 1/fb (2/fb for the Higgs 
search in cascade decays of sparticles) 198] 199 ^ . 

In contrast to this neutralino LSP scenario, the gravitino dark matter 
(GDM) scenario in the CMSSM is tightly constrained by the astrophysi- 
cal constraints on the cosmological abundances of light elements, as seen in 
Fig. H] [100J. However, such a scenario might have some advantages, e.g., 
by enabling the cosmological prediction for the abundance of 7 Li |1U1] to be 
improved, as also shown in Fig. Hfb). 

Recently, new attention has been focussed on the regions in which a 
metastable stau is the next-to-lightest sparticle (NSP) in a GDM scenario, 
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due to its ability to form bound states (primarily with 4 He). When such 
bound states occur, they catalyze certain nuclear reactions such as 4 He(D, 
7) 6 Li, which is normally highly suppressed due to the production of a low- 
energy 7, whereas the bound-state reaction is not |1U21 11031 1104] . In Fig. 
Ufa), the (mi/2,TOo) plane is displayed showing explicit element abundance 
contours jlOOj when the gravitino mass is m 3 / 2 = 0.2mo in the absence of stau 
bound-state effects. To the left of the solid black line the gravitino is not the 
LSP. The diagonal red dotted line corresponds to the boundary between a 
neutralino and stau NSP: above it, the neutralino is the NSP, and below it, 
the stau is the NSP. Very close to this boundary, there is a diagonal brown 
solid line. Above this line, the relic density of gravitinos from NSP decay is 
too high, i.e., 

^^n NSP h 2 > 0.12. (18) 
m NSP 

Thus we should restrict our attention to the area below this line. 

The very thick green line labelled 7 Li = 4.3 corresponds to the contour 
where 7 Li/H = 4.3 x 10 -10 , a value very close to the standard BBN result for 
7 Li/H. It forms a 'Vee' shape, whose right edge runs along the neutralino- 
stau NSP border. Below the Vee, the abundance of 7 Li is smaller than the 
standard BBN result. However, for relatively small values of m^, the 7 Li 
abundance does not differ very much from the standard BBN result: it is 
only when 7711/2 ^ 3000 GeV that 7 Li begins to drop significantly. The stau 
lifetime drops with increasing 777.1/2, and when r ~ 1000 s, at 77ii/ 2 ~ 4000 
GeV, the 7 Li abundance has been reduced to an observation-friendly value 
close to 2 x 10~ 10 as reported in [10511106] and shown by the (unlabeled) thin 
dashed (green) contours. 

The region where the 6 Li/ 7 Li ratio lies between 0.01 and 0.15 forms a 
band which moves from lower left to upper right. As one can see in the 
orange shading, there is a large region where the lithium isotopic ratio can 
be made acceptable. However, if we restrict to D/H < 4.0 x 10 -5 , we see 
that this ratio is interesting only when 7 Li is at or slightly below the standard 
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= 0.2m , tan (5 = 10 , |i > m 3/2 = 0.2m , tan (3 = 10 , |J. > 




m 1/2 (GeV) m 1/2 (GeV) 

Figure 4: The (m 1 / 2 ,iTi ) planes for m t = 172.7 GeV, A = 0, ji > 
and tan/3 = 10 with m 3 / 2 = 0.2m without (a) and with (b) the effects 
of metastable stau bound states included. The regions to the left of the solid 
black lines are not considered, since there the gravitino is not the LSP. In 
the orange (light) shaded regions, the differences between the calculated and 
observed light-element abundances are no greater than in standard BBN with- 
out late particle decays. In the pink (dark) shaded region in panel (b), the 
abundances lie within the ranges favoured by observation. The significances 
of the other lines and contours are explained in the text. 



BBN result. 

Turning now to Fig. H^b), we show the analogous results when the bound- 
state effects are included in the calculation. The abundance contours are 
identical to those in panel (a) above the diagonal dotted line, where the 
NSP is a neutralino and bound states do not form. We also note that the 
bound-state effects on D and 3 He are quite minimal, so that these element 
abundances are very similar to those in Fig. H^a). However, comparing panels 
(a) and (b), one sees dramatic bound-state effects on the lithium abundances. 
Everywhere to the left of the solid blue line labeled 0.15 is excluded. In 
the stau NSP region, this means that mi/2 <i 1500 GeV. Moreover, in the 
stau region to the right of the 6 Li/ 7 Li = 0.15 contour, the 7 Li abundance 
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drops below 9 x 10~ n (as shown by the thin green dotted curve). In this 
case, not only do the bound-state effects increase the 6 Li abundance when 
mi/a is small (i.e., at relatively long stau lifetimes), but they also decrease 
the 7 Li abundance when the lifetime of the stau is about 1500 s. Thus, 
at (m 1/2 ,m ) ~ (3200,400) GeV, we find that 6 Li/ 7 Li ~ 0.04, 7 Li/H ~ 
1.2 x 10" 10 , and D/H ~ 3.8 x 10" 5 . Indeed, when m 1/2 is between 3000- 
4000 GeV, the bound-state effects cut the 7 Li abundance roughly in half. 
In the darker (pink) region, the lithium abundances match the observational 
plateau values, with the properties 6 Li/ 7 Li > 0.01 and 0.9 x 10~ 10 < 7 Li/H < 
2.0 x 10~ 10 . This example demonstrates that it is possible to resolve the 
6 Li/ 7 Li by postulating GDM with a stau NSP. 

6 mSUGRA 

Minimal supergravity (mSUGRA) is often used as a basis for phenomeno- 
logical studies (301 E3l 1107] . The framework termed above the CMSSM is 
occasionally referred to as the mSUGRA. However, models based strictly on 
minimal supergravity should employ two additional constraints |108[ 1109] . 
One is a relation between the soft supersymmetry-breaking bilinear and tri- 
linear parameters: Bq = Aq — itlq, and the other is a relation between the 
gravitino and input scalar masses: 1713/2 = m®. In the simplest version of 
mSUGRA |110[ l30~t |23| 1107] . where supersymmetry is broken by a single 
field in a hidden sector, the universal trilinear soft supersymmetry-breaking 
terms are Aq = (3 — v3)too and bilinear soft supersymmetry-breaking term is 
Bq = (2 — \/3)mo, which is a special case of the general relation Bq = Aq— iiiq. 

Given such a relation between B and Aq, one can no longer use the 
standard CMSSM boundary conditions, in which mi/2, ttlq, Aq, tan/3, and 
sgn(fi) are input at the GUT scale, and then /1 and B are determined by 
the electroweak symmetry-breaking conditions. In this case, it is natural to 
use B as an input and calculate tan (3 from the minimization of the Higgs 
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potential [TtMfTUg] . 

Phenomenologically distinct planes may be determined by specifying a 
choice for Ao/mo, with the above-mentioned simplest hidden sector being 
one example. In Fig. two such planes are shown assuming m t = 172.7 
GeV with (a) A /m = 3 — \f?>, as predicted in the simplest model of su- 
persymmetry breaking jl 10] . and with (b) A /m = 2.0. We show in these 
(mi/2, mo) planes the contours of tan/3 as solid blue lines. Also shown are 
the contours where m x ± > 104 GeV (near-vertical black dashed lines) and 
m h > 114 GeV (diagonal red dash-dotted lines). The regions excluded by 
b — > S7 have medium (green) shading, those where the relic density of neu- 
tralinos lies within the WMAP range have light (turquoise) shading, and the 
region suggested by — 2 at 2-a has very light (yellow) shading, as in the 
CMSSM planes shown previously. As one can see, relatively low values of 
tan f3 are obtained in most of the visible planes. 

Another difference between the CMSSM and models based on mSUGRA 
concerns the mass of the gravitino. In the CMSSM, it is not specified and 
and can be taken suitably large so that the neutralino or the lighter stau is 
the LSP. In mSUGRA, the scalar masses at the GUT scale, mo, are deter- 
mined by (and equal to) the gravitino mass. In Fig. [5j the gravitino LSP 
and the neutralino LSP regions are separated by dark (chocolate) solid lines. 
Above these lines, the neutralino (or stau) is the LSP, whilst below them the 
gravitino is the LSP [111[ I112[ 1113] . As one can see by comparing the two 
panels, the potential for neutralino dark matter in mSUGRA models is de- 
pendent on Ao/mo- In panel (a), the only areas where the neutralino density 
are not too large occur where the Higgs mass is far too small or, at higher 
mo, the chargino mass is too small. At larger Ao/mo, the coannihilation strip 
rises above the neutralino-gravitino LSP boundary. In panel (b), we see the 
familiar coannihilation strip. It should be noted that the focus-point region 
is not realized in mSUGRA models as the value of \x does not decrease with 
increasing mo when Ao/mo is fixed and Bo = Aq — mo. There are also no 
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Figure 5: Examples of mSUGRA (m^jWo) planes with contours o/tan/3 
superposed, for > and (a) the simplest Polonyi model with Aq/tuq = 
3 — v^; and (b) A /m = 2.0, all with B = A — m . In each panel, we show 
the regions excluded by the LEP lower limits on MSSM particles and those 
ruled out by b — >■ sj decay (medium green shading): the regions favoured by 
g^ — 2 are very light (yellow) shaded, bordered by a thin (black) line. The dark 
(chocolate) solid lines separate the neutralino and gravitino LSP regions. The 
regions favoured by WMAP in the neutralino LSP case have light (turquoise) 
shading. The dashed (pink) line corresponds to the maximum relic density 
for the gravitino LSP, and regions allowed by BBN constraint neglecting the 
effects of bound states on NSP decay are light (yellow) shaded. 

funnel regions, as tan (3 is never sufficiently high. 

In the gravitino LSP regions, the NSP may be either the neutralino or 
stau, which are now unstable. (Note that in panel (b) there is also a re- 
gion which is excluded because the stau is the LSP.) The relic density of 
gravitinos is acceptably low only below the dashed (pink) line. This ex- 
cludes a supplementary domain of the (m^jMo) plane in panel (a) which 
has a neutralino NSP (the dotted (red) curve in panel (a) separates the 
neutralino and stau NSP regions). However, the strongest constraint is pro- 
vided by the effect of neutralino or stau decays on Big-Bang Nucleosynthe- 
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sis HH 1151 ESI HO ESI US]. Outside the light (yellow) shaded region, 
the decays spoil the success of BBN. 

7 Other Possibilities 

These cosmologically preferred regions move around in the (wii/2,mo) plane 
if one abandons the universality assumptions of the CMSSM. For example, if 
one allows the supersymmetry-breaking contributions to the Higgs masses to 
be non- universal (NUHM), the rapid-annihilation WMAP 'strip' can appear 
at different values of tan/3 and m^, as seen in Fig. [5] [ C2U\ 1121] . Rapid an- 
nihilation through the direct-channel H, A poles suppresses the relic density 
between the two parallel vertical WMAP strips at smaller values of m^, 
and the relic density is suppressed in the right-most strip because the neu- 
tralino LSP has a significant higgsino component. A complete exploration 
of the parameter space of the NUHM, which has two additional parameters 
compared to the CMSSM, lies beyond the scope of this review. 

The appearance of the (7711/2,7710) plane is also changed significantly if 
one assumes that the universality of soft super-symmetry-breaking masses 
in the CMSSM occurs not at the GUT scale, but at some lower renormal- 
ization scale [1221 11231 1124] , as occurs in some 'mirage unification' models 
[1251 11261 1127] . In this case, the sparticle masses are generally closer to- 
gether. As a consequence, the bulk, coannihilation, rapid-annihilation and 
focus-point regions approach each other and eventually merge as the mirage 
unification scale is reduced, as illustrated in Fig. [71 where they form an 'atoll'. 
At smaller values of the mirage unification scale, the atoll contracts and even- 
tually disappears, and there is no WMAP-compatible within the displayed 
portion of the (7711/2,7710) plane. In such 'GUTless' models, Vl^sph 2 falls be- 
low the WMAP range (1T7I) in larger regions of the (7711/2, 771q) plane than in 
the conventional CMSSM with unification at the GUT scale. 
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Figure 6: The (m^jMo) plane in the NUHM for tan (3 — 10, /i = 700 GeV 
and rriA = 400 GeV 1121^ . The colours of the shadings and contours are the 
same as in Fig. 



tan P = 10 , (X > , M in = 10 12 - 5 GeV 




m 1/2 (GeV) 



Figure 7: The (mi/2,mo) plane in a GUTless model with tan/3 = 10, fi > 
and v4 = 0, assuming universality at M in = 10 12 ' 5 GeV |123j . The colours 
of the shadings and contours are the same as in Fig. |2j 
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8 Summary 



As we have discussed above, there are many sound theoretical and phe- 
nomenological reasons to favour supersymmetric extensions of the Standard 
Model. In particular, supersymmetry predicts the existence of cold dark 
matter in a very natural way, and there are several plausible candidates for 
the lightest supersymmetric particle that would be present as a relic from 
the Big Bang. The most prominent candidate is the lightest neutralino, and 
we have described how its relic density may be calculated, and the regions of 
supersymmetric parameter space in which its density falls within the range 
favoured by astrophysics and cosmology. However, other candidates for the 
cold dark matter are also possible, such as the gravitino. In that case, the 
next-to-lightest supersymmetric particle would be metastable, and compar- 
isons between the observed light-element abundances and those predicted by 
Big-Bang Nucleosynthesis calculations impose important constraints on the 
parameter space. We have given examples of neutralino and gravitino dark 
matter scenarios in the minimal supersymmetric extension of the Standard 
Model, under various different theoretical assumptions. It will be for collider 
and dark matter detection experiments to determine which, if any, of these 
options has been adopted by Nature. 
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